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Studies of the phase diagram of the coupled sine circle map lattice have identified the presence of 
Q _ two dMinct .„™rs^.y classes of spatiote.por^ .tet.iUency vi. spat.o.e.po™! inter.it.ency 

^ ■ of the directed percolation class with a complete set of directed percolation exponents, and spatial 

" intermittency which does not belong to this class. We show that these two types of behavior 

are special cases of a spreading regime where each site can infect its neighbors permitting an 
initial disturbance to spread, and a non-spreading regime where no infection is possible , with 
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mapped on to an equivalent cellular automaton which shows a transition from a probabilistic 
' cellular automaton to a deterministic cellular automaton at the infection line. The origins of 



^0 , the spreading-non-spreading transition in the coupled map lattice, as well as the probabilistic to 

deterministic transition in the cellular automaton lie in a dynamical phenomenon, an attractor- 
^ ' widening crisis at the infection line. Indications of unstable dimension variability are seen in the 

neighborhood of the infection line. This may provide useful pointers to the spreading behavior 
seen in other extended systems. 
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I. INTRODUCTION 



The identification of tlie universality class of spatiotemporal intermittency Ij] in spatially 
extended systems has been a long standing problem in the literature. Early conjectures 
argued that the transition to spatiotemporal intermittency is a second order phase transition, 
and the transition falls in the same universality class as directed percolation j^. This 
conjecture has become the central issue in a long-standing debate , which is 

still not completely resolved. 

Studies of the coupled sine circle map lattice have thrown up a number of intriguing 
observations of relevance to this problem , [llj] . This system has regimes of spatiotem- 



poral intermittency (STI) with critical exponents which fall in the same universality class 
as directed percolation(DP), as well as regimes of spatial intermittency (SI) which do not 
belong to the DP class. Both these regimes lie on the bifurcation boundaries of the spa- 
tiotemporally fixed point solutions of the map. The spatiotemporally intermittent regime 
seen here has an absorbing laminar state, i.e. a laminar site remains laminar unless infected 
by a neighboring turbulent site. The burst states spread and can percolate through the 
entire lattice. The system shows a convincing set of DP exponents in this regime [o], Q, 
In the spatially intermittent regime, the laminar sites are frozen in time and the burst sites 
show temporally periodic or quasi-periodic behavior. The laminar sites do not get infected 
by neighboring turbulent sites. Hence, the spatially intermittent state is non-spreading and 
does not show DP exponents. Thus, both DP and non-DP behavior can be seen for different 
parameter regimes of the same system. 

In the present paper, we show that the infective DP behavior of spatiotemporal intermit- 
tency and the non-infective behavior of spatial intermittency are special cases of the more 
general spreading to non-spreading transition seen in this system. The spreading and non- 
spreading regimes are separated by a line which we call the infection line. Above the infection 
line, the burst states can infect neighboring laminar states and spread through the lattice, 
whereas below this line the burst states cannot infect their neighbors and the non-spreading 
regime is seen. The infection line intersects the bifurcation boundary of the synchronized 
solutions. Intermittent solutions are seen along this boundary, with the DP type of STI 
being seen above the infection line, and the non-DP SI being seen below the infection line. 
Spreading and non-spreading solutions are also seen off the bifurcation boundary. However, 
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the distribution of laminar lengths shows power-law scaling only for parameter values which 
are very close to the bifurcation boundary, and falls off exponentially as the parameter val- 
ues get more distant from the bifurcation boundary. Other exponents associated with DP 
behavior are also observed only along the bifurcation boundaries. 

Further insights into the spreading to non-spreading transition are obtained by map- 
ping the coupled map lattice (CML) onto a cellular automaton (CA). The spreading to 
non-spreading transition seen across the infection line maps on to a transition from a prob- 
abilistic cellular automaton to a deterministic cellular automaton. The dynamical origins 
of this transition lie in an attractor widening crisis which occurs at the infection line. The 
existence of this crisis can be inferred from the bifurcation diagram of the system. The 
finite-time Lyapunov exponents of the system fluctuate in the neighborhood of the infection 
line indicating the presence of unstable dimension variability. Thus the statistical charac- 
terizers of the system show signatures of a dynamical phenomenon. This result could have 
implications in a wider context. 

The paper is organized as follows. In section II, we give details of the coupled sine 
circle map lattice and the associated phase diagram and discuss the special behavior near 
the bifurcation boundary in these two regimes, namely spatiotemporal intermittency and 
spatial intermittency seen in the spreading and non-spreading regimes respectively. The 
signatures of the spreading and non-spreading regimes are also seen in the mapping of the 
coupled map lattice to a cellular automaton and its mean-field analysis. This is discussed 
in section III. In section IV, we discuss the dynamic origins of the two regimes, namely, the 
presence of an attractor-widening crisis and unstable dimension variability, as indicated by 
the fluctuating Lyapunov exponents at the infection line. We conclude with a discussion of 
the implications of our results. 



II. THE MODEL AND THE UNIVERSALITY CLASSES 



The coupled sine circle map lattice studied here is known to model the mode-locking 



behavior 



12l | seen in coupled oscillators, Josephson Junction arrays etc. The model is 



defined by the evolution equation 
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FIG. 1: shows the phase diagram of the coupled sine circle map lattice obtained with random initial 
conditions. The spreading and non-spreading regimes, separated by the infection line have been 
marked. Spatiotemporal intermittency of the directed percolation class is indicated by diamonds 
(O), whereas spatial intermittency has been marked with asterisks (*) and triangles (A). 

where i and t are the discrete site and time indices respectively and e is the strength of 
the coupling between the site i and its two nearest neighbors. The local on-site map, f{x) 
is the sine circle map defined as f{x) = x + — ^ sin(27rx), where, K is the strength of 
the nonlinearity and Q is the winding number of the single sine circle map in the absence of 
the nonlinearity. We study the system with periodic boundary conditions in the parameter 
regime < Q < ^ (where the single circle map has temporal period 1 solutions), < e < 1 
and K = 1.0. The phase diagram of this model evolved with random initial conditions is 
shown in Fig. [H 

Spatiotemporally fixed point solutions, in which all the sites relax to the fixed point 
X* = 2~: sin~^(^^), are seen in a large region of the phase diagram, and are indicated by 
dots. Two distinct regimes, separated by a line, which we call the infection line, can be 
identified in the phase diagram. A spreading regime is seen above the infection line, in 
which the burst states can infect their neighboring laminar states and spread through the 
lattice. The spreading nature of the bursts is strongly evident when a disturbance introduced 
in a laminar background eventually spreads to the entire lattice by infecting the neighboring 
laminar sites (Fig. [2](a)). This is in contrast to the non-spreading regime seen below the 
infection line, where the random initial conditions die down to bursts which are localized, 
and do not infect their neighboring laminar states (Fig. El^b)). 

In both these regions, spatiotemporal intermittency with co-existing laminar states ex- 
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FIG. 2: (Color online) shows the space-time plot of (a) spatiotemporal intermittency seen at = 
0.06, e = 0.7928 (one of the points marked with diamonds in Fig. [T|), and (b) spatial intermittency 
seen at = 0.047, e = 0.336 (marked with asterisk in Fig. [1]). The index i represents the lattice 
sites and the index t represents time. 

hibiting regular temporal behavior, and turbulent states with irregular temporal behavior, 
can be identified near the bifurcation boundary of spatiotemporally fixed point solutions in 
the phase diagram. The spatiotemporal intermittency seen in the spreading regime has been 
shown to belong to the directed percolation class. In the non-spreading regime, a different 
class of intermittency called the spatial intermittency is seen near the bifurcation boundary. 
We discuss these two classes of intermittency briefly in the following sections. 



A. Spatiotemporal intermittency of the DP class and the spreading regime 



Spatiotemporal intermittency belonging to the directed percolation class is seen in the 
vicinity of the bifurcation boundary of the spatiotemporally fixed point solutions in the 
spreading regime 0, Q, Some of the points that show this kind of spatiotemporal 

intermittency have been marked with diamonds (O) in the phase diagram (Fig. [1]). In 
this type of intermittency, the fixed point solution of the single sine circle map, x* acts 
as the laminar state, whereas the burst states lie in the [0, 1] interval. Since burst states 
cannot be spontaneously created, the laminar state acts as an absorbing state. The burst 
states can either 'percolate' through the lattice by infecting a neighboring laminar state, 
or die down to the laminar state. With time acting as the directed axis, the analogy with 
directed percolation is complete. Moreover, this class of spatiotemporal intermittency is 
remarkably free of coherent structures, some times called 'solitons', which spoilt the analogy 
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Static exponents 


Spreading Exponents 


z (3/vz (3 rj' C 


Tl 6 Ze 


STI-DP 


0.060 0.7928 
0.065 0.34949 
DP 


1.59(0.02) 0.17(0.02) 0.293 1.51(0.01) 1.68(0.01) 
1.59(0.03) 0.16(0.01) 0.273 1.50(0.01) 1.66(0.01) 
1.58 0.16 0.28 1.51 1.67 


0.315(0.007) 0.16(0.01) 1.26(0.01) 
0.303(0.001) 0.16(0.01) 1.27(0.01) 
0.313 0.16 1.26 


SI 


0.04 0.402 
0.047 0.336 


1.10(0.04) 
1.13(0.02) 





10| for definitions) obtained at 



TABLE I: The static and dynamic (spreading) exponents (see [9|, 
two of tlie directed percolation (DP) points (O-s) in Figure [1] are shown in the first two rows of 
the table. The universal DP exponents are listed in the third row. The laminar length distribution 
exponent, ^, calculated for spatial intermittency (SI) at two points marked by (A-s (quasi-periodic 
bursts) and (*-s (periodic bursts) respectively in Fig. [H are also listed. The error-bars are shown 
in the brackets. The data has been obtained for a lattice of size, N = 10^ and has been averaged 
over 10^ initial conditions. The laminar length distributions have been calculated for a lattice of 
size N = 10^ and averaged over 50 initial conditions. 



with directed percolation in other models such as the Chate Manneville coupled map lattice 

0. 

An entire set of static and dynamic scaling exponents was obtained in this parameter 
regime, which showed good agreement with the universal DP exponents. The exponents 
obtained, after averaging over 10^ initial conditions, at two such parameter values are listed 
in Table [B Similar exponents have been seen at points marked with diamonds {O) in the 
phase diagram (Figure [1]). These exponents also satisfy the hyperscaling relations for the 
static exponents (2/5/z/ = — 2 + r/'), and the spreading exponents (45 + 2rj = dzg) for 
d = 1. Hence, this class of spatiotemporal intermittency seen at the bifurcation boundary in 
the spreading regime belongs convincingly to the directed percolation class. These scaling 
exponents are however, seen only near the bifurcation boundary. The distribution of laminar 
lengths for the spreading solutions off the bifurcation boundary, show an exponential fall 
off. Thus STI of the DP class is a special case of the spatiotemporal behavior seen in the 
spreading regime. 

A different class of intermittency is seen below the infection line, i.e in the non-spreading 
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regime. This is the phenomenon of spatial intermittency. 



B. Spatial intermittency and the non-spreading regime 

Below the infection line, in the non-spreading regime, the laminar sites are the synchro- 
nized fixed point x* and the burst sites are either temporally frozen, periodic or aperiodic. 
These bursts are very stable and do not die down with time. 

In the vicinity of the bifurcation boundary, a special class of intermittency called spa- 
tial intermittency is seen. In this case, the temporal behavior of the burst states is either 
quasi-periodic or periodic. Some of these points have been marked with triangles (A, quasi- 
periodic) and asterisks (*, periodic) in the phase diagram (Fig. [1]). In this class of inter- 
mittency, the distribution of laminar lengths, obtained by averaging over different random 
initial conditions, scales as a power-law of the form P{1) ~ with ( as the associated 



scaling exponent [Uj. The exponent ( obtained in this case was found to be ~ 1.1 (Table 
[T]). A similar scaling exponent has been seen in the case of the spatial intermittency ob- 
served in inhomogeneous logistic map lattice |l^. The laminar length distributions seen for 
other non-spreading solutions at points off the bifurcation boundary show an exponential 
decrease. Thus the spatially intermittent solutions are special cases of the solutions in the 
non- spreading regime. 

Thus the sine circle coupled map lattice shows a transition from a spreading regime to 
a non-spreading regime at the infection line. In order to gain further insights into this 
transition, we map the coupled map lattice to a stochastic model, a probabilistic cellular 
automaton of the Domany-Kinzel type s], [l^ . 



III. MAPPING TO AN EQUIVALENT CELLULAR AUTOMATON 

Signatures of the spreading to non-spreading transition at the infection line are seen in 
a mapping of the coupled map lattice to a stochastic model namely, a probabilistic cellular 
automaton of the Domany-Kinzel type 0, |l^. The equivalent cellular automaton, set up to 
mimic the dynamics of the laminar and burst states in the coupled map lattice defined in 
eq. [H is defined on a one dimensional lattice of size A^. The state variable f * at site i and 
at time t is assigned the value f * = if the site is in the laminar state, and vj = 1 if the site 
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is in the burst state. 

As in the coupled map lattice equation [H the probability of the site i at time t + 1 
being in the burst state depends on the state of the sites i — l,i and i + 1 at time t. 
We therefore define the cellular automaton dynamics in this system by the conditional 
probability P{vl'^'^\vl_l,vl,vl^-^). There are 2^ possible configurations and the symmetry 
between the sites i — 1 and i + 1 in equation [1] gives the effective probabilities pkS as 
Po = P(1|000) , pi = P(l|001) = P(l|100), p2 = P(l|010), p3 = P(l|011) = P(l|110), 
P4 = P(l|101) and = P(l|lll). These probabilities then define the update rules of the 
cellular automaton. 

We estimate these probabilities pk, for a given set of parameter values, from the numerical 
evolution of the coupled map lattice with random initial conditions. The probabilities p^s 
are calculated by finding the fraction of sites i which exist in the burst state vj'^^ = 1 at 
time t + 1, given that the site i and its nearest neighbors i — 1 and i + 1 existed in state k 
at time t. That is, the probability pk is determined using pk = i^irrj^, where A^^q and 
are the number of sites, which at time t were the central sites of the configuration k, and at 
time t + 1 exist in the laminar states (f*"*""^ = 0) and the burst states (f*^""^ = 1) respectively. 
These probabilities, were extracted from a coupled map lattice of size = 10^ averaged 
over 14000 time steps discarding a transient of 1000 time steps, and averaged over 200 initial 





n 


e 


Po 


Pi 


P2 


P3 


P4 


P5 


S(DP) 


0.060 
0.073 


0.7928 
0.4664 


0.0 
0.0 


0.220 
0.150 


0.0 
0.0 


0.933 
0.938 


0.627 
0.439 


0.984 
0.993 


S 


0.070 


0.264 
0.248 


0.0 
0.0 


0.140 
0.050 


0.0 
0.0 


0.982 
0.989 


0.391 
0.160 


0.999 
0.999 


NS 


0.070 


0.232 
0.228 


0.0 
0.0 


0.000 
0.000 


0.0 
0.0 


1.000 
1.000 


0.000 
0.000 


1.000 
1.000 


NS(SI) 


0.031 
0.044 


0.420 
0.373 


0.0 
0.0 


0.000 
0.000 


0.0 
0.0 


1.000 
1.000 


0.000 
0.000 


1.000 
1.000 



TABLE II: shows the probabilities p^s obtained in the spreading (S), non-spreading (NS) regimes 
and at directed percolation (DP) and spatial intermittency (SI) points. 
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FIG. 3: (Color online) shows the probabilities associated with the update rules , , p4 , plotted 
as a function of the coupling strength, e at = 0.065. In the non-spreading regime (e < Cic), the 
probabilities are either or 1, whereas they lie in the interval (0, 1) in the spreading regime (e > Cic). 



conditions. The probabilities pk obtained in the spreading and non-spreading regimes in the 
phase diagram are listed in Table [Tll 

It is clear that the condition for an absorbing state is realized in the probability po = 
P(1|000) which is seen to be zero in both regimes. The probabilities pi = P(l|001) = 
P(l|100) and p4 = P(l|101) essentially define infection probabilities, which estimate the 
probability of a laminar site being infected by its burst neighbor or neighbors, to change to 
a burst site. 

We can see from Table [III that these probabilities show drastically different behavior in 
the spreading and non-spreading regimes. In the case of the spreading regime, the prob- 
abilities obtained are seen to lie in the open interval (0,1). Therefore, the dynamics in 
the spreading regime is described by a probabilistic cellular automaton (PCA) wherein the 
cellular automaton rules are probabilistic in nature. In contrast, in the case of the non- 



spreading regime, in addition to po and p2 which are zero for the STI of the DP type |15| . 
the infection probabilities pi and p^ are also seen to be equal to zero. Hence, the infection 
probabilities characterize the non-infective nature of the bursts. Moreover, the probabilities 
P3 = P(l|011) = P(l|110) and ps = P(l|lll) take the value 1, which indicates the robust- 
ness of the burst states. We note that the probabilities obtained here only take the values or 
1. Hence, we obtain a deterministic cellular automaton (DCA) in the non-spreading regime, 
wherein given a state k at time t, the state of the site i at time t + 1 is decisively known with 
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probability zero or one. This is further illustrated in Figure [3l where the probabilities have 
been plotted as a function of the coupling strength e at Q = 0.065. The probabilities show a 
distinct change to values other than zero or one, at the point where spreading regime starts 
viz. e > (See inset in Fig. 1). Thus a transition from a probabilistic to a deterministic 
cellular automaton is seen at the infection line. 

It is therefore clear that the spatiotemporal intermittency of the directed percolation 
class seen in the spreading regime can be modeled by a PCA, whereas spatial intermittency 
observed in the non-spreading regime can be modeled by a DCA. In the following section, we 
show that the probabilistic cellular automata obtained at the STI of the DP points indeed 
mimics the dynamics of the coupled map lattice and shows scaling exponents matching with 
the universal directed percolation exponents. 

A. Probabilistic cellular automaton and directed percolation 

As mentioned previously, the dynamics seen in the case of spatiotemporal intermittency 
in the spreading regime resembles that of directed percolation. The stochastic nature of the 
burst states is further emphasized in the probabilistic cellular automaton obtained in this 
regime. In this section, we show that the probabilistic cellular automata obtained at points 
in the phase diagram which show DP-like spatiotemporal intermittency, indeed show scaling 
behavior which match with the directed percolation class. 

Figure H] shows some of the physical quantities associated with the static and dynamical 
properties of the system. The scaling exponents obtained for these quantities are also shown 
in the figure. The probabilities associated with the PCA obtained at the parameter value 
Q = 0.06, e = 0.7928, were found to be po = 0.0, pi = 0.213116, p2 = 0.000081,^3 = 
0.902515,^4 = 0.569464 and = 0.967595. When the cellular automaton is prepared with 
random initial conditions, such that the sites take either or 1 values with equal probabilities 
at time t = 0, the burst state (1 state) is seen to die down with time. The relaxation time r, 
which determines the time taken by the lattice to relax to a laminar state (0 state), scales as 
a function of the size of the lattice L as r ~ (Fig. Hl^a)). The order parameter m, defined 
as the fraction of burst sites in the lattice, decreases with time as m ~ 
In contrast, when a pair of burst states are introduced in a lattice prepared in the laminar 
state at time t = 0, the burst states are seen to grow with time. The dynamical quantities 
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FIG. 4: (Color online) shows the log-log (base 10) plot of physical quantities calculated for the 
probabilistic cellular automaton associated with the DP point 0, = 0.06, e = 0.7928. The figure 
shows the (a) escape time r vs length of the lattice L, (b) order parameter m vs time t, and the 
spreading properties (c) the fraction of burst sites N(t) vs t, and (d) the survival probability P{t) 
plotted as a function of time t. The fits to the plots and the corresponding scaling exponents (error 
bars in brackets) obtained are also shown. 

associated with the growth of the burst states viz. the fraction of burst sites N{t), and the 
survival probability P{t), which is defined as the fraction of initial conditions which yield a 
non-zero number of burst states at time t, are shown in Fig. W^c) and (d). The respective 
scaling exponents obtained show good agreement with the directed percolation class. Thus 
the cellular automaton model obtained in the region of spatiotemporal intermittency seen 
in the spreading regime exhibits scaling behavior similar to the directed percolation class, 
as expected. 

We note that in the case of the non-spreading regime, the deterministic cellular automaton 
obtained acts like an identity mapping. Therefore, the scaling behavior in the laminar length 
distribution depends strongly on the initial condition namely, the fraction of burst states 
in the lattice. To mimic the dynamics exhibited by the coupled map lattice, we calculate 
the fraction of burst sites in a coupled map lattice of size 20000 aX Q = 0.04, e = 0.402, 



11 




FIG. 5: shows the laminar length distribution obtained for the deterministic cellular automaton 
obtained at the spatial intermittency point 0, = 0.04, e = 0.402. The exponent obtained is C = 
1.106 ±0.005. 

after discarding 60000 transients and averaging over 50 random initial conditions. The 
fraction of burst sites came out to be p = 0.0684. Fig. [5] shows the scaling behavior of the 
distribution of laminar lengths for the deterministic cellular automaton, obtained for this 
initial condition, in which bursts were introduced with probability p. The exponent seen in 
the cellular automaton turns out to have the value ~ 1.1 as in the coupled map lattice. 

A mean-field study of the cellular automaton gives further insights into the transition 
from probabilistic to deterministic cellular automaton, as well as into the initial condition 
dependence of the cellular automata. We discuss this in the subsequent section. 



B. Mean-field analysis of the cellular automaton 

We study a mean field approximation of the probabilistic cellular automaton which gives 
a better understanding of the probabilistic to deterministic cellular automaton transition 

n 

as well as gives probability bounds for the directed percolation-like behavior [16|]. Let rrit 
and m^+i be the density of burst states in the lattice at the i:*'* and t + 1*'* time step. The 
mean-field for the probabilistic cellular automaton can be defined, according to the cellular 
automaton rules, as 

mt+i = {2pi + P2)mt{l - mtf + (2j93 + Pi)ml{l - rrit) + p^ml (2) 
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FIG. 6: (Color online) shows the stability regions of the three fixed points of the mean-field 
equation O m^,m2, and rrtg, in the {^1,^3,^4} phase space. The probabilities of the DP points 
(marked with diamonds), as well as for points in the spreading regime at = 0.065 (dotted line) 
have been shown. The triangle corresponds to the probabilities in the non-spreading regime. 

By approximating = I and p2 = 0, the mean field equation reduces to 

mi+i = 2pi mt{l - mtf + (2j93 + p^) mf{l - rrit) + (3) 

The three fixed points of this equation for a set of parameters {^1,^3,^4} are m\ = 
0,m2 = i-|-2p^f-2p\-p4 ' "^3 ~ ^- ^ linear stability analysis of the mean-field equation 
indicates that the fixed point m\ = 0, which corresponds to the absorbing state, is stable 
when pi < |. The fixed point is stable when 2p^ + P4 < 2 and 2pi > 1. The fixed point 
mg = 1, which corresponds to the completely turbulent state, is stable when 2p^ + P4> 2. 
The stability regions of these three fixed points in the {pi,P3,P4} space have been indicated 
in Fig. El The 2p3 + p^ = 2 plane, above which = 1 solutions are found to be stable, is 
shown in the figure. 

We see that both = and = 1 solutions are stable, when pi < | and 2p3 -|-p4 > 2. 
This is the co-existence region. The system settles down to either the completely laminar 
state {m\ = 0) or the completely turbulent state (rrig = 1) in the co-existence region, 
depending on the initial condition rrit aX t = 0. The solutions settle down to the m\ = 
solution if initial density, mo G [0,7712) whereas, they tend towards mg = 1 if tuq G (m2, 1]. 
For instance, we iterate the map equation [3] by choosing p4 = 0.9 and study the {m,pi,p2,} 
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FIG. 7: shows the {?7i,Pi,P3} phase diagram of the mean field map at p4 = 0.9 with initial m at 
time t = chosen as (a) niQ = 0.1 and (b) rriQ = 0.9. In the co-existence region (2p3 > 1.1 and 
Pi < 0.5), the density of burst sites m, settles down to either m\ or rrig depending on the initial 
density mo chosen. 

phase diagram for two initial conditions mo = 0.1 and mo = 0.9. The {m,pi,p-^} phase 
diagram obtained is shown in Figure [7](a) and (b) respectively. We see that the solutions 
settled down to m* = in the co-existence region in the first case, whereas they settle down 
to m* = 1 in the latter case. Hence, the solutions in the co-existence region are initial 
condition dependent. 

Interestingly, we find that the probabilities associated with the spreading regime of the 
coupled map lattice, including those obtained for the directed percolation points in the 
spreading regime, lie in the co-existence region (Figure E]). Due to this, the probabilistic 
cellular automaton has a very strong initial condition dependence in the spreading region. 
The correct choice of initial conditions leads to the laminar absorbing state. Other choices 
can end up easily in the m = 1 state. However, the points at which directed percolation 
is seen are exceptions to this. The probabilistic cellular automata obtained at these points 
show scaling behavior consistent with the directed percolation phenomena irrespective of 
the choice of initial conditions, as was discussed in the previous section. 

The probabilities associated with the deterministic cellular automaton seen in the non- 
spreading regime (i.e. pi = 0, ps = 1, P4 = 0) lie at the vertex belonging to the co-existence 
region in this cube. This point has been marked with a triangle in Figure [61 As we approach 
the spatial intermittency points along the bifurcation boundary, the cellular automaton 
probabilities obtained at the directed percolation points tend towards this vertex. Similarly, 
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the probabilities obtained in the spreading regime tend towards this vertex, as we cross 
the infection hne and enter the non-spreading regime. Therefore, the deterministic cellular 
automaton is a limiting case of the probabilistic cellular automaton seen in the spreading 
regime. 

This confirms that the presence of both spreading and non-spreading regimes separated 
by the infection line in the phase diagram, which accounts for the presence of both directed 
percolation-like spatiotemporal intermittency and spatial intermittency in the phase dia- 
gram, is manifested in the form of a transition from probabilistic to deterministic cellular 
automaton at the infection line. However, the origins of the spreading nature of the burst 
states above the infection line are not clear. In the next section, we identify the dynamical 
origins of the spreading states. 



IV. CRISIS AND UNSTABLE DIMENSION VARIABILITY AT THE INFECTION 
LINE 



It is clear from the preceding discussion that the behavior of the burst states changes 
drastically at the infection line leading to the spreading to non-spreading transition and the 
existence of distinct classes of spatiotemporal intermittency in the phase diagram. We now 
have indications that the origin of the spreading states lies in an attractor-widening crisis 
at the infection line as well as the presence of unstable dimension variability in the vicinity 
of the infection line, as we will see in this section. 

A crisis is said to have occurred in a dynamical system when a sudden change in the 
attractor takes place when a system parameter is changed We observe such a phe- 

nomenon in our coupled map lattice, when we change the strength of the coupling between 
sites, e, for a given Q in the vicinity of the infection line. This has been illustrated in Fig. 
m Here, the variable x associated with a typical site has been plotted over 500 time steps, 
as a function of the coupling strength e at Q = 0.065. 

The range of e values on the vertical axis of Fig. [1] cuts across the infection line at 
etc = 0.252 for Q = 0.065. The bifurcation diagram clearly shows that an attractor widening 



crisis 



ISl] appears at this point. Similar behavior is seen for other sites. The spreading regime 



seen in the phase diagram emerges exactly at the point at which the attractor widens, with 
the non-spreading regime corresponding to the pre-widening regime. This widening also 
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FIG. 8: shows the bifurcation diagram of the coupled map lattice in which the variable x associated 
with a typical site has been plotted over 500 time steps, as a function of the coupling strength e in 
the neighborhood of the infection line, at = 0.065. 

identifies the point at which the equivalent cellular automaton undergoes a probabilistic 
to deterministic transition. In the pre-crisis region, each site follows either a periodic or 
quasi-periodic trajectory and is not infected by the behavior of its neighbors. Thus, its CA 
analogue is deterministic as listed in Table [Tll In the post-crisis regime, each site is able 
to access the full x range, as well as infect its neighbors, and the bursting and spreading 
behavior characteristic of the spreading regime is seen. This is reflected in the equivalent 
cellular automaton by a transition to probabilistic behavior (Table HTl) . It is to be noted 
that the volume of the attractor in phase space will be much larger post-crisis, as compared 
to the pre-crisis volume. 

We also see indications of unstable dimension variability in the vicinity of the infection 
line. An attractor is said to have unstable dimension variabilit y, i f it possesses periodic 



orbits with different number of stable and unstable directions |19l. l20[|. A trajectory visiting 
the neighborhood of these periodic orbits, experiences a fluctuating number of unstable 
directions, as time progresses. A signature of this phenomenon can be seen in the finite time 
Lyapunov exponents (FTLE) of the system which fluctuate about zero, as the dynamics 
evolves. 

For the coupled map lattice of eq. [U, which is an N-dimensional map of the form x — > f (x) , 
let the Jacobian evaluated at the initial condition xq be denoted as Df(xo). The time-n 
Lyapunov exponents (FTLE) for this system, are obtained by diagonalizing the Jacobian 
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FIG. 9: (Color online) shows (a) the fraction of positive finite-time Lyapunov exponents (time=20) 
as a function of time t obtained at = 0.065, e = 0.254, and (b) the distribution of the largest 
finite-time Lyapunov exponent (time= 100) plotted for 0, = 0.065 and e = 0.245,0.252,0.260. 

Df"(xo) of the n-times iterated map f", the k^'^ time-n Lyapunov exponent being defined 

as 

Afc(xo,n) = ^logAfc (4) 

where, is the fc*^ eigenvalue of the time-n Jacobian Df"(xo). 

We find that the finite-time Lyapunov exponents obtained for our system exhibit fluctu- 
ations about zero in the vicinity of the infection hne. This can be seen in Fig. M^a), where 
the fraction of positive time-20 Lyapunov exponents obtained at = 0.065, e = 0.254, is 
plotted as a function of time. The fraction of positive exponents is seen to vary with time, 
indicating that a fraction of the Lyapunov exponents fluctuate about zero, as the dynam- 
ics evolves. Hence, the trajectory of the coupled map visits periodic orbits with different 
number of unstable directions, in the vicinity of the infection line. 

This is further illustrated in Fig. MJo), which shows the distribution of the largest time- 
100 Lyapunov exponent, obtained a.t Q = 0.065 for various values of e. The distribution 
has been calculated by collecting data for a lattice of size N = 200 over 10000 time steps 
and by averaging over 200 initial conditions. We find that the distribution of the largest 
time-100 Lyapunov exponent is constrained to the negative values for parameters lying in 
the non-spreading regime (eg. e = 0.245 in Fig. [9]^b)). For parameters close to the infection 
line (viz. e = 0.252), the distribution shifts towards the positive side, until it lies entirely in 
the positive side of the axis for parameters above the infection line (eg. e = 0.260). Hence, 
the distribution of Lyapunov exponents obtained at e = 0.252 illustrates further that the 



17 



largest finite-time Lyapunov exponent vacillates around zero near the infection line. 

Therefore, we infer from Fig. [9](a) and (b), that the finite-time Lyapunov exponents 
fluctuate about zero, near the infection line, thereby confirming the presence of unstable 
dimension variability at the infection line. Hence, we see that the spreading nature of the 
bursts above the infection line can be attributed to the attractor-widening crisis as the well 
as the presence of unstable dimension variability at the infection line. 



V. DISCUSSION 



To conclude, the spatiotemporal intermittency of the directed percolation class and the 
spatial intermittency of the non-directed percolation class, seen along the bifurcation bound- 
aries of spatiotemporally fixed point solutions, are special cases of the spreading and non- 
spreading regimes seen off the bifurcation boundaries. The two regimes are separated by 
the infection line which intersects with the bifurcation boundary of the fixed point solutions 
at the point where the cross-over between the directed percolation and non-directed per- 
colation behavior takes place. Thus the behavior seen in coupled sine circle map lattice is 
organized around the locations of the bifurcation boundaries of the fixed point solutions, 
and the infection line. The existence of two distinct universality classes, in the phase di- 
agram of the sine circle map lattice is refiected in the transition of the equivalent cellular 
automaton from the probabilistic phase to the deterministic phase and the concomitant sup- 
pression of the spreading or infectious modes. The dynamic origins of this transition lie in 
an attractor-widening crisis. We also find evidence for the existence of unstable dimension 
variability in the neighborhood of the infection line. Thus, in this system, there is a direct 
connection between a dynamical phenomenon viz. a crisis in an extended system and the 
statistical properties of the extended system viz. the exponents and universality classes. 
Similar directed percolation to non-directed percolation transitions have been seen in other 
coupled map lattices, as well as in pair contact processes, solid on solid models and models of 
non-equilibrium wetting 2l(]. Our results may have useful pointers for the analysis of other 
systems, and thus contribute to the on-going debate on the identification of the universality 
classes of spatiotemporal systems. 
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